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Project  Summary 


This  project  focuses  on  the  benefits  and  challenges  associated  with  the  utilization  of  Earth-Moon  (lunar) 
libration  point  orbits  (LPOs)  for  communications  with  and  surveillance  from  military  satellites.  DoD- 
relevant  topical  areas  include  Space  Situational  Awareness  and  Operationally  Responsive  Space 
capabilities  with  emphasis  on  Space-Based  Space  Surveillance  (SBSS),  orbital  camouflage,  and  reducing 
risks  of  both  purposeful  interference  and  collisions  with  space  debris  that  are  associated  with  space  assets 
in  low  Earth  (LEO)  or  geosynchronous  (GEO)  orbit.  In  the  first  year  of  the  project  efforts  concentrated  on 
1)  finding  fuel-optimal  invariant  manifold-based  transfers  via  impulsive  maneuvers  to  the  lunar  Li  and  L3 
Lagrange  points,  2)  studying  Earth  and  GEO  observability  from  single  and  multiple  LPOs,  and  3)  analysis  of 
thrust  dispersion  errors  at  manifold  injection  and  the  development  of  guidance  laws  to  mitigate  their 
effects.  In  the  second  year  efforts  concentrated  on  1)  finding  fuel-optimal  transfers  via  impulsive 
maneuvers  to  a  selection  of  halo  orbits  about  the  lunar  Li  and  L3  LPOs,  2)  finding  fuel-optimal  transfers  to 
a  Liapunov  orbit  about  the  lunar  Li  point  using  mixed  impulsive  and  continuous  thrust,  3)  studying  Earth 
and  GEO  observability  from  single  LPO  about  the  Li  and  L3  points  and  a  constellation  of  spacecraft  in 
multiple  Li  and  L3  LPOs,  and  4)  the  development  of  continuous  (LQR-based)  and  impulsive  guidance  laws 
to  mitigate  the  effects  of  thrust  errors  for  the  optimal  impulsive  manifold-based  transfers  to  the  Li  LPO. 
In  the  third  year  efforts  have  concentrated  on  1)  obtaining  improved  impulsive  transfers  from  Earth  to  a 
halo  orbit  about  the  lunar  L3  point,  2)  obtaining  statistical  fuel  budgets  for  impulsive  guidance  for  use  in 
fuel-optimal  transfers  to  lunar  Li  LPOs,  3)  obtaining  low-cost  periodic-gain  feedback  control  laws  for 
station  keeping  about  lunar  Li  LPOs  and  comparing  their  performance  with  more  traditional  control 
schemes.  In  total  ten  conference  papers  [1-10]  and  one  published  journal  paper  [11]  resulted  from  this 
investigation,  as  well  as  three  more  currently  in  review  or  preparation  [12-14].  In  addition,  portions  of 
the  results  from  this  project  were  presented  in  two  talks  [15-16]  by  the  PI,  one  of  which  was  given  at  the 
Air  Force  Research  Laboratories,  Kirtland  AFB,  Albuquerque.  The  following  subsections  outline  the  main 
contributions. 

Technical  Background 

A.  Problem  Formulation 

The  Circular  Restricted  Three-Body  Problem  (CR3BP)  models  the  motion  of  a  particle  of  negligible  mass 
under  the  influence  of  two  larger  bodies,  or  primaries  (Figure  1).  The  larger  of  the  bodies  is  termed  the 
primary  and  the  smaller  is  termed  the  secondary.  The  primaries  rotate  in  circular  orbits  about  the  system 
barycenter.  Additionally,  the  reference  frame  rotates  about  the  barycenter  at  the  same  rotation  rate  as 
the  two  primaries.  The  x-axis  extends  from  the  origin  through  the  secondary,  the  z-axis  extends  in  the 
direction  of  the  angular  momentum  of  the  system,  and  the  y-axis  completes  the  right-hand  coordinate 
frame.  The  reader  is  directed  to  Szebehely  [17]  for  a  derivation  of  the  equations  of  motion.  The  Jacobi 
Constant,  C,  is  the  only  integral  of  motion  which  emerges  in  the  rotating  system. 

B.  Lagrange  Points 

The  CR3BP  has  five  equilibrium  points  with  coordinates  (xe,ye,ze)  called  Lagrange  or  libration  points  and 
denoted  by  Li-L5.  These  points  can  be  found  by  setting  the  derivatives  in  the  CR3BP  equations  equal  to 
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zero.  Using  this  procedure,  one  can  see  that  for  three  of  them,  called  the  collinear  libration  points  L1-L3, 
the  ye  and  ze  components  are  zero.  For  the  Earth-Moon  system  with  gravitational  constant  p=0. 01215, 
the  liberation  points  are  located  from  the  barycenter  with  the  following  dimensionless  coordinates  (for 
distances  in  km  multiply  by  the  average  Earth-Moon  distance  384,400  km):  Li:  (xe,ye,ze)=(0.8369,0,0), 
L2:(l. 156,0,0),  and  L3:(-1.005,0,0).  Two  other  equilibrium  points  L4  and  L5  (called  triangular  points)  are 
located  at  (xe, ye, ze)=(0.48785, ±0.86603,0)  .  All  five  Lagrange  points  of  Earth-Moon  system  have  been 
shown  in  Figure  1. 

C.  Periodic  Orbits  about  Lagrange  Points 

Families  of  periodic  and  quasi-periodic  orbits  which  exist  about  the  libration  points  have  been  studied 
extensively  by  many  researchers  [18].  The  quasi-periodic  orbits  around  liberation  points  are  called 
Lissajous  orbits  which  wind  around  a  torus  but  never  repeat  [19].  Halo  orbits  are  a  special  case  of  Lissajous 
orbits  where  the  in-plane  and  out-of-plane  frequencies  are  equal  thus  producing  periodic  motion,  and 
were  discovered  by  Farquhar  [24].  Large  amplitude  lunar  halo  orbits  appear  to  form  a  halo  about  the 
Moon  when  viewed  from  the  Earth  [20,  25].  Halo  orbits  are  three-dimensional;  a  two-dimensional 
periodic  orbit  in  the  moon's  orbital  plane  is  referred  to  as  a  Lyapunov  orbit  [21-22]. 

Several  methods  exist  for  numerically  computing  libration  point  orbits  (LPOs).  The  orbits  presented  here 
were  computed  using  a  Single-Shooting  Algorithm  [23].  The  focus  of  this  project  is  on  using  Halo  orbits  at 
the  Li  and  L3  points.  Figure  2  shows  views  of  Halo  orbits  around  Li  in  3D  and  four  different  planes.  Four 
views  of  northern  Earth-Moon  L3  halo  orbits  are  shown  in  Figure  3.  Orbits  at  L3  can  be  totally  planar  (e.g., 
Lyapunov  orbits),  or  have  a  very  large  vertical  component.  Additionally,  certain  orbits  at  L3  actually  have 
perigees  that  are  less  that  the  geosynchronous  radius  (42,184  km). 

D.  Stability  of  Libration  Points  and  Halo  Orbits 

To  determine  the  stability  of  the  libration  points,  a  perturbed  displacement  is  defined  and  substituted  into 
the  equations  of  motion  which  are  then  Taylor-expanded  up  to  linear  terms  about  (xe,ye,ze).  Evaluating 
the  characteristic  equations  for  the  collinear  Lagrange  points  (Li_3)  reveals  their  saddle-type  instability 
while  L4-5  maintain  linear  stability.  For  the  Li  libration  point  in  the  Earth-Moon  system,  the  eigenvalues  are 
A  =  ±2.9321,  ±2.3344i,  ±2.2688i  which  correspond  to  the  Es,  Eu,  and  Ec  subspaces,  respectively.  For  the  L3 
libration  point  in  the  Earth-Moon  system,  the  eigenvalues  are  A  =  ±0.7548,  ±1.0038,  ±1.8914i.  To  obtain 
the  stability  of  halo  orbits,  Floquet  theory  must  be  used.  The  monodromy  matrix  is  defined  as  the  state 
transition  matrix  computed  at  the  orbital  period,  and  its  eigenvalues  (Floquet  multipliers)  can  be  used  to 
obtain  a  stability  index.  A  stable  orbit  has  a  stability  index  of  one,  and  stability  indices  greater  than  one 
corresponds  to  unstable  orbits.  As  the  stability  index  increases,  the  orbital  stability  decreases.  Halo  orbits 
at  Li  and  L2  can  have  stability  indices  that  vary  greatly,  depending  on  their  size  and  proximity  to  the 
secondary.  For  example,  in  the  Earth-Moon  system,  a  family  of  simple  Li  halo  orbits  can  have  stability 
indices  that  range  from  1  to  over  1000.  However,  halo  orbits  at  L3  are  relatively  much  more  stable,  with 
stability  indices  ranging  from  1  to  1.382.  For  periodic  orbits  in  the  CRTBP,  it  is  useful  to  introduce  the 
parameter  t  ,  which  is  used  to  describe  the  position  of  a  particle  on  a  periodic  orbit.  This  is  similar  to  the 
mean  anomaly  in  the  two-body  problem,  and  is  therefore  a  measure  of  the  time  elapsed  on  the  periodic 
orbit  and  not  an  angular  measurement.  The  value  of  t  ranges  from  0-360  degrees,  increasing  in  the 
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direction  of  orbital  motion,  t  is  defined  to  be  zero  at  the  orbit's  initial  x-axis  crossing  in  the  +y  direction, 
and  -180  at  the  x-axis  crossing  in  the  -y  direction,  as  shown  in  Figure  4. 

E.  Invariant  Manifolds 

Stable  (unstable)  manifolds  consist  of  the  set  of  all  possible  trajectories  a  spacecraft  could  follow  as  it 
asymptotically  approaches  (departs)  a  libration  point.  The  nonlinear  stable  (Ws)  (unstable  (Wu))  manifolds 
can  be  found  by  numerically  integrating  the  CR3BP  equations  with  an  initial  condition  of  the  state  of  each 
libration  point  orbit  perturbed  in  the  direction  of  its  stable  (unstable)  eigenvector  backward  (forward)  in 
time.  The  size  of  perturbation  applied  to  produce  the  stable  (unstable)  manifolds  in  Earth-Moon  system 
is  10A-4  in  non-dimensional  units.  The  stable  and  unstable  manifolds  of  a  Li  halo  orbit  can  be  seen  in 
Figure  5  as  tube-like  structures.  Each  orbit  has  two  associated  stable  and  unstable  manifold  sets:  one 
corresponding  to  a  positive  perturbation,  and  one  corresponding  to  a  negative  perturbation.  The  stable 
and  unstable  manifolds  of  periodic  orbits  can  be  found  at  any  point  along  the  orbit.  To  find  the  manifolds 
at  each  point  along  the  orbit,  it  is  first  necessary  to  propagate  the  stable  and  unstable  eigenvectors  to  the 
desired  position  along  the  orbit  using  the  state  transition  matrix.  The  unstable  (stable)  manifold  for  a 
periodic  orbit  can  be  calculated  by  propagating  a  state  perturbed  in  the  unstable  (stable)  direction  forward 
(backward)  in  time.  The  stability  of  an  orbit  dictates  how  fast  a  particle  can  travel  through  the  manifolds. 
Orbits  with  a  lower  stability  index  (more  stable  orbits)  have  stable  (unstable)  manifolds  which  takes  longer 
time  for  a  particle  to  get  close  (far)  to  the  periodic  orbit.  As  a  comparison,  stable  manifolds  of  two  LPOs 
around  Li,  one  with  v=2.7  and  the  other  with  v=1031.3,  have  been  propagated  for  75  days  in  Figure  6.  The 
stable  manifold  of  the  more  unstable  LPO  gets  close  the  vicinity  of  the  Earth,  but  the  manifold 
corresponding  to  the  one  with  lower  stability  index  is  still  in  the  vicinity  of  LPO.  Halo  orbits  around  L3  are 
much  more  stable  then  the  halo  orbits  around  Li  (the  stability  index  for  L3  Halo  orbits  are  less  than  1.4); 
therefore,  the  related  manifolds  for  L3  halo  orbits  are  much  slower  and  it  takes  very  long  times  to  get  close 
to  L3  halo  orbits  using  these  trajectories.  Figure  7  shows  that  for  stable  manifolds  of  L3  halo  orbits,  after 
one  year  propagation  there  are  few  perigees  less  than  20,000  km  and  after  two  year  propagation  there 
are  some  perigees  less  than  20,000  km  which  all  occur  after  540  days  of  propagation.  As  a  result  of  this 
study,  it  was  found  that  invariant  manifolds  are  not  practical  to  design  trajectories  for  transfers  to  L3  halo 
orbits. 

F.  Pseudo-Manifolds 

In  order  to  take  advantage  of  the  convergent/divergent  behavior  of  manifolds  and  construct  transfer 
trajectories  from  Earth  which  will  arrive  to  the  orbit  in  a  reasonable  time  frame,  the  concept  of  a  Pseudo- 
Manifold  was  introduced  in  the  first  year  of  the  project.  The  Pseudo-Manifold  varies  from  the  traditional 
manifold  in  two  distinct  ways: 

1)  The  perturbation  to  construct  the  trajectory  is  only  given  to  the  velocity  components  of  the  state,  and 

2)  the  magnitude  of  the  perturbation  is  much  larger  (1000  times  larger  than  the  perturbation  needed  to 
produce  the  true  manifolds).  Figure  8  illustrates  the  difference  between  the  stable  manifold  and  the 
pseudo-stable  manifold  for  a  L3  halo  orbit  with  Jacobi  Constant  C  =  2.4207.  The  location  on  the  orbit  to 
give  the  perturbation  for  producing  the  manifold  trajectories  was  selected  as  t=  17.  The  applied 
perturbation  (e)  for  the  stable  manifold  is  10A-4  in  non-dimensional  units,  and  for  the  pseudo-stable 
manifold,  £  =  130  m/s  =0.1269  non-dimensional  units.  The  stable  manifold  trajectory  was  propagated  for 
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312  days  while  the  pseudo-stable  manifold  trajectory  was  propagated  for  only  48  days.  The  pseudo-stable 
manifold,  however,  quickly  departs  the  vicinity  of  the  halo  orbit  and  performs  flybys  of  both  the  Earth  and 
the  Moon.  Over  the  given  time  frames,  the  perigee  of  the  stable  manifold  is  74,454  km  and  occurs  on  day 
293  of  the  trajectory  while  the  pseudo-stable  manifold  has  a  perigee  of  7,056  km  on  day  46  of  the 
trajectory. 

G.  Earth  Coverage  and  GEO  Surveillance  from  Libration  Point  Orbits 

Another  goal  of  this  project  is  to  study  the  visibility  of  the  Earth's  surface  as  well  as  the  geosynchronous 
(GEO)  belt  using  a  constellation  of  spacecrafts  in  halo  orbits  about  the  Li  and  L3  libration  points.  Space- 
Based  Space  Surveillance  (SBSS)  is  essential  to  improve  and  maintain  our  knowledge  of  the  many  objects 
in  space  such  as  debris  and  other  spacecraft.  Unlike  the  LEO  and  GEO  regions,  there  is  no  other  spacecraft 
in  Earth-Moon  halo  orbits  and  consequently  there  is  no  danger  of  collision  with  orbital  debris  and  reduced 
risk  of  intentional  interference.  Also,  SBSS  is  not  affected  by  weather  or  the  time  of  day,  so  it  can  record 
useful  data  more  consistently  than  ground-based  systems.  Furthermore,  because  observations  made  from 
space  do  not  have  an  atmosphere  diffusing  what  little  light  is  available,  SBSS  has  the  potential  to  see 
dimmer  objects.  Past  SBSS  missions  have  used  telescopes  in  sun-synchronous  orbits  [15].  SBSS  has  been 
proven  useful  by  many  studies  and  multiple  delivered  spacecraft.  Some  notable  examples  are  the  U.S.  Air 
Force's  Space-Based  Visible  (SBV)  sensor  on  the  MSX  spacecraft,  the  Space-Based  Space  Surveillance 
(SBSS)  Block  10  System  built  by  Boeing  and  Ball  Aerospace  [26],  and  the  Canadian  Sapphire  satellite  [27]. 
The  SBV  sensor  was  part  of  the  first  system  to  perform  SBSS  [28]. 

Summary  of  the  Most  Important  Results 

A.  Optimal  Transfers  to  Li,  L3,  and  Li  Halo  Orbits 

Two  optimization  methods  were  used  for  finding  minimum-fuel  manifold-based  co-planar  transfers  to  Li 
and  L3:  a  grid  search  and  a  genetic  algorithm.  The  grid  search  optimization  method  is  a  direct  search  type 
algorithm  that  can  be  computationally  inexpensive  when  the  number  of  variables  in  the  cost  function  is 
minimal  [29].  In  this  method,  a  grid  is  constructed  in  the  decision  space  and  the  cost  function  is  evaluated 
at  each  grid  point.  The  decision  space  includes  all  possible  values  of  each  variable.  The  grid  point  that 
minimizes  the  AVtotai  is  therefore  the  optimal  solution.  The  variable  0  is  assumed  to  vary  from  0  degrees 
to  359  degrees  by  increments  of  1  degree.  The  resulting  optimal  trajectory  for  transfer  from  a  300  km  LEO 
to  Li  occurs  at  0=3°  with  AVtotai=3.72  km/s  and  a  time  of  flight  of  24.8  days,  while  the  optimal  trajectory 
to  the  L3  point  takes  155.9  days  with  AVtotaF3.49  km/s  and  0=3°. 

For  finding  optimal  two-impulse  transfer  trajectories  to  a  LI  halo  orbit,  a  specific  southern  halo  orbit  with 
a  Jacobi  constant  of  C=3.10,  stability  index  of  216.06,  period  of  12.1  days,  and  z-amplitudes  of  29,714  km 
and  -42,592  km  is  selected  and  an  initial  300  km  LEO  or  35,780  km  GEO  is  assumed.  This  initial  orbit  is 
considered  to  be  inclined  relative  to  the  Moon's  (x-y)  orbital  plane  with  a  free  inclination  that  is  solved  for 
in  the  optimization  process  in  the  case  of  departing  LEO  or  constrained  as  20  deg  from  the  x-y  plane  in  the 
case  of  departing  GEO,  where  0  is  the  departure  angle  measured  counter-clockwise  from  the  x-axis.  The 
initial  AV  for  a  transfer  trajectory  that  intersects  the  stable  manifold  of  Li  or  L3  along  with  the  angle  0  are 
the  variables  need  to  be  determined  through  the  optimization.  A  second  impulsive  maneuver  is  needed 
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to  change  the  velocity  of  the  spacecraft  so  that  it  matches  the  velocity  of  the  stable  manifold  at  the 
intersection  point.  A  minimum  distance  function  is  used  to  find  the  closest  point  for  which  the  transfer 
and  manifold  intersect,  along  with  the  corresponding  velocities.  By  altering  the  values  of  each  variable, 
an  infinite  number  of  transfers  can  be  produced  to  simulate  the  transfer  of  a  spacecraft  from  LEO  to  the 
Li  or  L3  stable  manifold.  A  combination  of  a  genetic  algorithm  (GA),  primer  vector  theory,  and  a  level  I 
differential  corrector  were  used  to  solve  for  the  free  variables  that  produce  an  optimal  two-impulse  low- 
cost  transfer,  where  the  population  is  described  by  four  variables:  the  point  on  the  halo  orbit  from  which 
the  stable  manifold  is  propagated  (x),  the  time  the  manifold  is  propagated  for  (ts),  the  pseudo  right 
ascension  to  the  ascending  node  (pseudo-Q,  defined  relative  to  the  x-y  plane),  and  the  total  delta-v,  which 
is  the  dependent  variable.  The  target  position  on  the  stable  manifold  is  the  final  point  after  which  it  has 
been  propagated  for  time  ts.  Primer  vector  theory  is  then  used  in  the  fitness  function  to  determine  the 
optimal  departure  point  from  LEO  or  GEO  by  allowing  for  an  initial  coast  that  improves  the  cost  function. 
The  GA  continues  for  100  generations  or  until  there  has  been  no  improvement  in  the  dependent  variable 
for  ten  generations. 

The  optimal  trajectories  obtained  by  this  method  are  shown  in  Figs.  9-10,  in  which  the  GA  scatter 
plots  are  shown  converging  to  the  optimal  solution.  In  Fig.  9  for  LEO  departure  the  AVtotai  for  this  transfer 
is  3.61  km/s  with  a  AVi  =  3.07  km/s  and  AV2  =  0.54  km/s.  This  transfer  is  described  by  its  x  value  of  198 
deg,  pseudo-Q  of  155.07  deg,  and  a  ts  of  34.68  days.  The  algorithm  found  an  optimal  initial  LEO  inclined 
21.82  deg  relative  to  the  Moon's  orbital  plane  while  the  time  of  flight  is  37.93  days.  In  Fig.  10  for  GEO 
departure  the  AVtotai  for  this  transfer  is  1.14  km/s  with  a  AVi  =  0.89  km/s  and  AV2  =  0.25  km/s.  This  transfer 
is  described  by  its  x  value  of  280  deg,  pseudo-Q  of  176.03  deg,  and  a  ts  of  35.64  days.  The  time  of  flight 
is  40.75  days.  By  comparison,  direct  transfers  and  alternative  locally  optimal  transfers  from  LEO  and  GEO 
are  shown  in  Figs.  11-12,  respectively.  The  direct  transfer  from  LEO  requires  a  AVtotai  of  3.77  km/s,  which 
is  160  m/s  more  than  that  for  the  optimal  transfer.  However,  the  direct  transfer  takes  only  3.82  days.  On 
the  other  hand,  a  locally  optimal  transfer  that  requires  less  travel  time  and  a  small  amount  of  extra  delta- 
v  is  shown  where  AVtotai=3.63  km/s  and  the  transfer  time  is  29.56  days.  The  direct  transfer  from  GEO 
requires  a  AVtotai  of  1.44  km/s,  which  is  300  m/s  more  than  that  for  the  optimal  transfer.  However,  the 
direct  transfer  takes  only  5.12  days.  On  the  other  hand,  a  locally  optimal  transfer  that  requires  less  travel 
time  and  a  small  amount  of  extra  delta-v  is  shown  that  requires  10  fewer  days  than  the  globally  optimal 
transfer  and  requires  50  m/s  more  in  AVtotai. 

B.  Analysis  of  Thrust  Dispersion  Effects 

Thrust  dispersion  analyses  for  the  manifold  injection  maneuvers  (assuming  impulsive  burns)  were 
examined  for  the  Li  halo  orbit  optimal  transfer  discussed  above.  Monte  Carlo  analyses  of  1000  simulations 
each  were  performed  with  effects  of  thrust  magnitude  and  direction  errors  analyzed  separately  and  in 
combination  on  the  resulting  miss  distance  and  transfer  times.  Gaussian  dispersions  for  the  same  three 
errors  are  introduced  into  the  system:  one  for  the  AV2  magnitude  (assuming  la  for  1%  error)  and  two  for 
the  AV2  direction  in  the  burn-plane  and  out-of-plane  directions  (assuming  la  for  1  degree  error).  The 
results  shown  in  Figs.  13-14  for  departure  from  LEO  and  GEO,  respectively,  illustrate  the  need  for  a 
guidance  scheme  to  mitigate  the  effects  of  thrust  errors.  For  the  combined  thrust  magnitude  and 
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direction  errors  with  LEO  departure,  the  mean  miss  distance  was  17678.06  km  with  a  la  of  3196.48  km 
while  the  transfer  times  had  a  mean  of  22.22  days  with  a  la  of  1.66  days.  For  the  combined  thrust 
magnitude  and  direction  errors  with  GEO  departure,  the  mean  miss  distance  was  6884.71  km  with  a  la 
of  2583.56  km  while  the  transfer  times  had  a  mean  of  23.50  days  with  a  la  of  0.77  days.  Note  the  presence 
of  outliers  in  the  miss  distance  dispersions  representing  trajectories  that  completely  miss  the  halo  orbit 
by  thousands  of  kilometers.  These  large  final  position  errors  at  the  halo  orbit  insertion  point  do  not  allow 
a  ballistic  transfer  as  desired. 

C.  LQR-Based  Continuous  Thrust  Guidance 

A  neighboring  optimal  guidance  scheme  based  on  the  linear  quadratic  regulator  (LQR)  was  implemented 
assuming  the  spacecraft  has  a  mechanism  for  low  continuous  thrust  during  its  travel  time  on  the  stable 
manifold  in  order  to  correct  for  the  thrust  dispersion  errors  and  reduce  the  miss  distance  at  the  halo  orbit 
insertion  point.  The  effectiveness  of  the  guidance  control  for  the  transfers  from  LEO  and  GEO  is  seen  in 
Figs.  15-16  which  shows  the  transfer  plot,  position  and  velocity  state  differences  (as  compared  with  the 
nominal  states  on  the  stable  manifold),  and  the  three  elements  of  the  control  input  vector.  For  the  LEO 
departure  case  with  guidance  in  Figure  15,  the  guided  position  and  velocity  deviations  which  are  now 
significantly  reduced.  The  final  miss  distance  is  3.05  m,  a  significant  improvement  over  the  cases  of 
transfer  without  guidance  and  the  optimal  transfer  without  any  injection  errors.  This  shows  that  small 
targeting  errors  that  lead  to  significant  miss  distances  can  be  compensated  for  with  the  guidance  scheme. 
The  guidance  controller  uses  a  total  delta-v  of  15.365  m/s  during  its  34.57  day  trip  on  the  manifold.  It  can 
be  seen  that  there  is  large  control  magnitude  immediately  after  the  delta-v2  maneuver  which  quickly 
decays  for  the  rest  of  the  transfer.  The  maximum  output  thrust  used  in  this  control  scheme  for  this  transfer 
is  100.0067  m/sA2.  This  thrust  output  is  feasible  with  the  currently  existing  low-thrust  engines.  For  the 
GEO  departure  case  with  guidance  in  Figure  16,  the  final  miss  distance  is  1.96  m.  The  guidance  controller 
uses  a  total  delta-v  of  7.052  m/s  during  its  35.64  day  trip  on  the  manifold.  The  maximum  output  thrust 
used  in  this  control  scheme  for  this  transfer  is  0.003042  m/sA2.  Again,  this  output  thrust  is  feasible  with 
currently  existing  low-thrust  engines.  Finally,  we  repeat  the  GEO  departure  case  while  limiting  the  thrust 
magnitude  used  in  the  guidance  scheme  to  a  maximum  of  Tm  =  0.01  mm/sA2.  Figure  17  shows  the 
corresponding  results  which  show  that  the  state  deviations  from  those  on  the  manifold  are  still  reduced 
to  zero  but  over  a  longer  period  than  in  Fig.  16.  The  final  miss  distance  in  this  case  is  2.5  m  while  the 
controller  uses  a  total  delta-v  of  12.795  m/s  during  its  35.64  day  trip  on  the  manifold.  It  is  noticed  that 
the  cost  is  increased  from  the  case  without  thrust  limits. 

D.  Impulsive  Guidance 

A  more  realistic  guidance  strategy  is  to  use  impulsive  trajectory  correction  maneuvers  (TCMs)  instead  of 
the  continuous  thrust  that  results  from  the  LQR-based  solution.  Our  impulsive  guidance  strategy  utilizes 
the  state  transition  matrix  (STM)  to  propagate  state  deviations.  Two  impulsive  TCMs  (constituting  a  "burn 
pair")  are  used:  the  first  burn  at  time  ti  re-targets  the  manifold  and  the  second  burn  at  time  t2  corrects 
the  velocity  error  upon  arrival  at  the  manifold.  Control  strategies  include  options  where  only  single  burns 
or  full  burn  pairs  are  implemented.  The  timing  of  the  TCMs  was  obtained  first  by  a  brute  force  trial-and- 
error  approach  and  then  by  using  the  Liapunov  Exponents  (LEs)  associated  with  the  halo  orbit  stable 
manifold.  The  LEs  associated  with  the  stable  manifold  portion  of  the  optimal  two-impulse  transfers  in 
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Figs.  9-10  for  departure  from  LEO  and  GEO  are  shown  in  Figure  18,  in  which  the  maximum  values 
correspond  to  perigees  and  the  minimum  values  correspond  to  apogees.  The  initial  hypothesis  that  the 
TCMs  should  be  performed  at  apogee  locations  (minimum  values  of  LEs)  was  confirmed  through  testing 
multiple  ideas  on  how  to  use  them  in  implementing  the  TCM  burn  times.  The  resulting  transfers  from  LEO 
and  GEO  with  single  impulsive  TCMs  calculated  from  the  LEs  are  shown  in  Figs.  19-20,  in  which  the  thrust 
errors  (1%  magnitude  and  1  deg  in-plane  and  out-of-plane)  are  incorporated  in  the  manifold  injection 
maneuver  as  well  as  in  all  subsequent  TCMs,  which  are  indicated  by  the  black  circles  in  the  transfer  plot 
with  the  relative  magnitudes  shown  in  the  plots  below.  In  the  case  of  the  transfer  from  LEO,  the  miss 
distance  at  halo  orbit  insertion  is  reduced  to  4.858  km  while  the  transfer  time  is  34.569  days.  The  delta- 
v  cost  for  the  impulsive  guidance  is  22.034  m/s.  In  the  case  of  the  transfer  from  GEO,  the  miss  distance 
at  halo  orbit  insertion  is  reduced  to  1.092  km  while  the  transfer  time  is  35.642  days.  The  delta-v  cost  for 
the  impulsive  guidance  is  8.634  m/s.  This  shows  that  the  impulsive  guidance  scheme  is  effective  in 
mitigating  the  position  errors  on  the  stable  manifold  as  well  as  the  halo  insertion  miss  distance  that  result 
from  thrust  errors  in  the  manifold  injection  maneuver  and  the  TCMs  themselves  compared  to  the  case 
where  no  guidance  is  used. 

E.  Optimal  Mixed  Impulsive/Continuous  Thrust  Transfers  to  Li  Liapunov  Orbits 

Optimal  transfer  trajectories  were  designed  for  a  spacecraft  using  mixed  impulsive  and  continuous  thrust 
propulsion  to  depart  low-Earth  orbit  and  enter  a  specified  planar  Lyapunov  orbit  at  the  interior  LI  Earth- 
Moon  Lagrange  point  in  the  framework  of  the  planar  circular  restricted  three-body  problem.  Due  to  the 
fact  that  impulsive  and  continuous  thrust  transfers  have  their  own  relative  advantages  and  disadvantages, 
the  combination  of  the  two  propulsion  types  has  been  suggested  as  a  way  to  combine  th  benefits  of  both 
systems  into  transfer  trajectory  designs.  The  flight  time  and  impulsive/continuous  thrust  weighting  factor 
are  specified  in  advance,  while  the  results  from  two  separate  performance  indices  based  on  minimum  fuel 
and  minimum  control  energy  were  compared.  The  continuous  dynamic  optimization  problem  was 
reformulated  as  a  discrete  optimization  through  direct  transcription  and  collocation,  which  then  was 
solved  using  nonlinear  programming  software.  The  optimal  transfer  trajectory  results  were  analyzed  using 
the  Pareto  front  solutions.  Fig.  21  shows  direct  and  spiral  departure  transfer  trajectories  and  the 
corresponding  control  acceleration  histories  using  a  quadratic  performance  index  in  the  optimal  control 
problem  (OCP).  The  results  for  varying  numbers  of  nodes  are  shown  to  illustrate  when  convergence  is 
achieved:  for  the  direct  departure  case  (which  includes  a  mandatory  coast  period  after  the  initial 
impulsive  maneuver  to  depart  GEO)  197  nodes  is  shown  to  be  sufficient  while  for  the  spiral  departure 
case  99  nodes  is  sufficient  (for  which  the  control  history  is  smooth).  By  contrast,  Fig.  22  shows  direct  and 
spiral  departure  transfer  trajectories  and  the  corresponding  control  acceleration  histories  using  a 
minimum  fuel  performance  index  for  spiral  departure  in  the  OCP,  in  which  99  nodes  is  shown  to  be 
sufficient  for  convergence  and  the  thrust  history  consists  of  two  periods  of  maximum  thrust  at  the  pre¬ 
set  thrust  limit  (dashed  red  line)  and  two  periods  of  coast  arcs  with  no  thrust.  The  bang-off-bang  thrust 
profile  was  found  to  be  explained  in  the  context  of  primer  vector  theory  when  the  performance  index 
minimizes  the  actual  fuel  used  as  opposed  to  the  case  of  a  quadratic  performance  index. 

F.  Efficient  Transfers  to  L3  Halo  Orbits 

Work  has  continued  on  finding  efficient  transfers  from  LEO  to  halo  orbits  about  the  Earth-Moon  L3 
Lagrange  point.  Due  to  very  slow  nature  of  the  L3  halo  orbit's  stable  manifolds,  Davis  et  al.  [5]  proposed 
transfers  based  on  L3  halo  orbit  pseudo-manifolds  (defined  by  applying  a  large  perturbation  in  the  velocity 
component  of  the  eigenvector)  from  LEO  to  lunar  L3  halo  orbits  to  take  advantage  of  the 
convergent/divergent  behavior  of  manifolds  and  construct  transfer  trajectories  from  Earth  which  will 
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arrive  to  the  orbit  in  a  reasonable  time  frame.  For  transfer  from  a  185  km  LEO  to  a  northern  L3  halo  orbit 
with  C  =  2.42  and  period  of  27.2  days,  Fig.  23  shows  a  pseudo-manifold-based  transfer  trajectory  found  in 
this  way  where  the  insertion  point  at  the  halo  orbit  is  given  by  t  =  17.36°  and  the  velocity  perturbation 
(also  the  amount  of  delta-v2  for  the  2-impulse  transfer)  is  £  =  130  m/s.  The  first  maneuver  is  a  tangential 
maneuver  of  3.07  km/s  and  is  executed  in  LEO  to  inject  onto  the  pseudo-manifold.  The  total  AV  for  this 
transfer  is  3.2  km/s  while  the  total  transfer  time  is  45.8  days.  In  general  the  total  delta-v  cost  is  between 
3.17  to  3.28  km/s  with  flight  times  between  43  and  84  days.  However,  the  lowest  delta-v2  that  was 
reported  is  not  less  than  130  m/s. 

In  such  mission  designs  with  two  impulsive  burns,  it  is  important  that,  for  a  given  total  AV,  as  much  as 
possible  should  be  used  for  the  first  burn  departing  LEO  while  the  second  burn  should  be  as  small  as 
possible.  This  way,  the  amount  of  fuel  the  spacecraft  should  carry  to  the  halo  orbit  is  minimized.  For  this 
purpose,  recent  efforts  have  investigated  alternative  strategies  that  result  in  more  desirable  transfers. 
Using  a  grid  search  applied  to  360  different  points  of  t  (every  one  degree)  where  at  each  r  location  on  the 
halo  orbit,  the  orbit's  nominal  velocity  is  perturbed  in  32400  different  directions,  departure  from  a  variety 
of  LEO  altitudes  is  achievable.  Figure  24  shows  the  achievable  LEO  altitudes  that  result  from  a  100  m/s 
perturbation  to  the  nominal  state  on  an  L3  halo  orbit.  The  locations  for  the  perturbations  are  denoted  by 
theiri  value.  The  halo  orbit  has  a  Jacobi  constant  of  C=2.3981  and  period  of  27.13  days.  The  transfer  time 
was  limited  to  110  days.  Transfers  were  constructed  for  each  of  the  trajectories  with  a  LEO  altitude  less 
than  400  km.  It  was  assumed  that  the  first  maneuver  to  depart  from  LEO  is  a  tangential  maneuver  onto 
the  trajectory  in  order  to  minimize  the  value  of  AVi.  The  total  AVs  of  the  resulting  transfer  range  between 
2.94  and  3.14  km/s.  The  total  AV  increases  as  the  LEO  altitude  increases.  The  trajectories  found  in  the 
Figure  24  can  have  a  time  of  flight  between  45  to  110  days.  Figure  25  shows  different  families  of  transfers 
from  a  185  km  altitude  LEO  to  the  same  L3  LPO  in  which  the  AV2,  AVtotai,  and  time  of  flight  (TOF)  are 
indicated  for  each  family.  A  common  feature  shared  between  all  families  shown  in  Figure  25  is  the 
inclusion  of  a  lunar  flyby  with  an  altitude  of  a  few  thousand  kilometers,  while  a  feature  that  distinguishes 
between  different  families  is  the  number  of  the  Earth  flybys  which  is  a  measure  for  the  time  of  flight 
ranges  between  the  families  of  the  transfers.  The  number  of  the  Earth  flybys  varies  between  no  flybys  to 
three  flybys  if  the  time  of  flight  is  limited  to  120  days.  Another  interesting  observation  one  can  find  in 
Figure  25  is  that  different  families  are  distinguishable  even  between  categories  which  have  similar  ranges 
of  TOF.  As  an  example,  families  of  transfers  in  Figures  b,  f,  i  all  have  one  Earth  flyby  but  their  TOF  varies 
between  three  different  ranges:  40  days  to  47  days  for  the  first  family,  62  to  65  for  the  families  in  the 
Figure  f,  and  62  to  65  days  for  the  families  of  the  transfers  showed  in  Figure  i.  The  same  feature  is 
observable  for  the  families  which  perform  no  Earth  flyby  and  those  who  have  two  Earth  flybys.  Figure  26 
shows  the  C3  energy  respect  to  the  halo  orbit  insertion  maneuver  for  the  constructed  transfers  using 
colored  dots.  Different  transfers  colored  from  blue  to  red  as  their  transfer  time  increases.  C3  varies 
between  -1.9  and  -2.18.  There  is  no  correlation  between  C3  and  AV2and  also  between  C3  and  time  of  flight, 
while  a  distinguishable  trend  can  be  found  for  AV2  and  time  of  flight;  as  AV2  is  decreased  the  time  of  flight 
increases. 

The  three  parameters:  AV2,  C3  energy,  and  ToF,  all  effects  on  the  cost  and  feasibility  of  the  transfers. 
C3  effects  for  the  launch  vehicle  and  AV2  effects  on  the  type  of  engine  spacecraft  need  for  the  second  burn 
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and  also  the  amount  of  fuel  it  needs  to  carry.  Therefore,  it  is  necessary  to  find  a  trade-off  between  these 
three  quantities.  A  transfer  requiring  the  smallest  insertion  maneuver  into  the  halo  orbit  can  be 
constructed.  Figure  27  shows  such  a  transfer  which  the  minimum  possible  AV2=15  m/s  for  insertion  into 
the  halo  orbit.  This  transfer  requires  an  escape  energy  of  C3=-2.05  kmA2/sA2  and  takes  119.6  days  to  insert 
into  the  halo  orbit  while  performing  a  lunar  flyby  and  only  one  Earth  flyby.  On  the  other  hand  a  transfer 
can  be  designed  to  have  the  shortest  time  of  flight.  Figure  28  shows  this  transfer  which  takes  only  31  days 
for  the  spacecraft  to  insert  into  the  L3  halo  orbit.  The  C3  escape  energy  for  this  transfer  is  -2.02  kmA2/sA2 
while  it  needs  90  m/s  burn  to  insert  into  the  halo  orbit. 

G.  GEO  Observability  Analysis  from  Libration  Point  Orbits 

Two  steps  were  defined:  Visibility  analysis  of  the  GEO  belt  using  Matlab  programs,  and  validation  of  the 
results  using  NASA's  Ephemeris  model  in  the  orbit  analysis  software  Satellite  Took  Kit  (STK).  As  the  first 
step,  a  study  was  conducted  to  find  the  visibility  coverage  using  several  constellations  of  spacecrafts  in  Li 
and  l_3  Halo  orbits  of  Earth-Moon  system.  By  analyzing  various  constellations,  this  study  tried  to  find  a 
balance  between  coverage  and  cost  to  provide  coverage  for  largest  percentage  of  the  Earth  surface  while 
minimizing  the  number  of  spacecrafts  in  the  proposed  constellation.  There  are  also  some  other  factors 
most  be  considered  such  as  orbit  insertion  costs,  orbit  maintenance  costs,  and  link  budget  considerations. 
For  this  purpose,  a  grid  was  defined  by  latitude/longitude  points  for  equal  dispersion  along  a  sphere.  Then, 
at  each  point  along  the  halo  orbit,  using  the  line-of-sight  criteria,  it  was  determined  which  points  of  this 
grid  are  visible  to  spacecraft  in  halo  orbits,  and  finally  four  percentages  were  calculated  as  a  parameters 
of  visibility.  These  parameters  are  defined  as:  Total  percentage  Visible:  The  total  percent  of  the  Earth's 
surface  visible  during  one  full  halo  orbital  period  (assuming  no  rotation  of  the  Earth);  Maximum 
percentage  Visible:  The  percent  of  the  Earth's  surface  visible  at  the  optimal  viewing  location  on  the  orbit; 
Minimum  percentage  Visible:  The  percent  of  the  Earth's  surface  visible  at  the  worst  viewing  location  on 
the  orbit;  Average  percentage  Visible:  The  average  percentage  of  the  Earth's  surface  visible  during  an 
orbital  period.  To  validate  the  results  of  the  first  step,  the  analyses  used  a  combination  of  STK  and 
software  written  in  MATLAB.  Libration  point  orbits  were  generated  in  MATLAB  and  saved  as  STK 
ephemeris  files.  The  LPO  ephemeris  files  were  then  imported  into  STK  using  the  stkExternal  propagator. 
Other  satellite  orbits,  such  as  sun-synchronous  orbits  and  GEO  orbits,  were  imported  into  STK  using  the 
available  catalog  of  two-line  elements  (TLEs).  Finally,  STK's  J4Perturbation  integrator  was  used  to  model 
any  arbitrary  orbit  that  was  used  apart  from  catalog  orbits. 

A  grid  of  1502  points  can  be  defined  by  latitude/longitude  points  for  equal  dispersion  along  sphere 
with  N=4°  in  latitude  and  M=7.2°  in  longitude.  Then,  at  each  point  along  the  halo  orbit,  using  the  line-of- 
sight  criteria,  it  should  be  determined  which  points  of  this  grid  are  visible  to  spacecraft  in  halo  orbit  and 
finally  a  percentage  were  calculated  as  a  value  of  visibility.  The  orbit  is  propagated  for  four  orbital  periods 
of  L3  orbit  (approximately  108  days).  2400  different  orbit  combinations  are  tested  in  this  study.  With  one 
spacecraft  on  E-M  Li  halo  orbit,  approximately  50%  of  the  Earth's  surface  at  a  given  time  can  be  observed. 
Most  of  the  E-M  Li  halo  orbits  provide  similar  coverage  characteristics  and  If  the  Earth  were  fixed  (i.e.,  no 
rotation),  less  than  60%  of  the  Earth's  surface  could  be  seen  in  1  full  halo  period;  therefore,  Halo  orbits  at 
Li  can  only  see  the  "near  side"  of  the  Earth.  Figure  29  shows  the  four  percentage  parameter  for  a 
spacecraft  on  Li  halo  orbit.  Also,  a  spacecraft  on  L3  halo  orbit,  on  average,  can  see  over  97%  of  the 
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geosynchronous  region  and  there  is  only  small  blockage  due  to  Earth.  Figure  30  illustrates  the  maximum 
amount  of  geosynchronous  region  visible  during  1  orbital  period  in  green  color  (best  case)  and  the 
minimum  amount  of  the  geosynchronous  region  visible  during  1  orbital  period  (i.e.,  the  portion  blocked 
by  the  Earth)  as  the  worst  case.  Different  constellation  like  four  different  combinations  of  two  spacecrafts 
in  two  different  orbits  (1  s/c  in  L3  orbit  (Northern  or  Southern)  and  1  s/c  in  Li  orbit  (Northern  or  Southern)) 
was  examined  and  finally  with  a  constellation  of  two  spacecraft  in  L3  halo  orbit  (one  in  a  northern  and 
another  in  a  southern  one)  and  one  spacecraft  in  a  northern  Li  halo  orbit,  there  can  be  a  coverage  of  up 
to  99.8  percent  of  Earth  surface  (as  seen  in  Figure  51)  and  all  examined  constellations  can  see  at  least  97% 
of  Earth  surface  on  average. 

To  generate  the  halo  orbit  in  STK,  first,  an  approximate  orbit  was  generated  using  a  single-shooting 
algorithm  in  the  Circular  Restricted  Three-Body  Problem.  The  approximate  orbit  was  converted  into  a 
continuous,  ballistic  trajectory  in  the  high-fidelity  model  of  the  solar  system  using  a  multiple-shooting 
differential  corrector  [18-19].  The  high-fidelity  model  uses  the  gravity  of  the  Sun,  Moon,  and  all  planets, 
modeled  as  point-masses  traveling  in  orbits  approximated  by  the  Jet  Propulsion  Laboratory's  Planetary 
and  Lunar  Ephemerides  DE405  [20].  No  assumptions  were  made  regarding  solar  radiation  pressure  due 
to  the  early  analyses  of  these  trajectories,  but  the  trajectories  will  be  adjusted  in  future  analyses  to 
compensate  for  solar  radiation  pressure.  Figure  32  illustrates  the  Li,  L2,  and  L3  trajectories  used  in  the 
analyses  presented  here.  The  brute  force  analysis  was  performed  for  an  observing  spacecraft  orbiting 
about  L3.  Looking  at  the  whole  set  of  possible  orbits,  we  see  some  interesting  patterns  appear  in  Figure 
33.  At  low  altitudes,  the  target  spacecraft  is  visible  for  an  extremely  small  percentage  of  the  total  time. 
Above  5,000  km,  the  viewing  opportunities  become  much  more  abundant,  but  there  are  only  a  few  orbits 
of  interest  above  that  point -namely,  GPS  and  GEO.  At  high  altitudes,  the  limiting  factor  for  target  visibility 
is  the  phase  angle  -  the  viewer  only  sees  the  illuminated  side  of  the  target  about  22.2%  of  the  time.  This 
follows  from  the  assumption  made  that  the  phase  angle  must  be  less  than  40°  for  visibility  (80°/360°  = 
22.2%).  The  range  of  altitudes  that  houses  the  GPS,  GLONASS,  and  Galileo  navigation  constellations  is 
sparsely  populated,  so  GEO  emerges  as  the  best  orbital  regime  for  observation  from  a  libration-point 
orbit.  Figure  33,  Figure  34,  and  Figure  35  summarize  these  findings.  It  was  observed  that  for  low-altitude 
orbits,  the  visibility  was  severely  affected  by  the  Earth  being  either  directly  in  front  of  or  directly  behind 
the  target  satellite.  LEO  spacecraft  also  spend  much  more  time  in  eclipse  than  higher  orbits,  resulting  in 
the  limited  time  visible  seen  above  in  Figure  32.  The  effect  of  each  constraint  as  a  function  of  altitude  is 
shown  below  in  Figure  36. 

To  examine  the  visibility  of  GEO  from  a  LPO,  a  complete  analysis  was  run  independently  for  observer 
spacecraft  at  Li,  L2,  and  L3.  Figures  37,  38,  and  39  show  the  percentage  of  GEO  spacecraft  that  are  visible 
from  an  Li,  L2,  and  L3  orbit  over  time,  respectively.  Figure  40  shows  the  percentage  of  the  total  number  of 
accesses  that  are  shorter  than  a  given  duration.  Figure  41  is  a  histogram  showing  the  frequency  of  various 
access  durations.  There  are  two  main  spikes  in  the  number  of  access  durations  for  Li:  one  at  approximately 
2.25  hours  and  one  between  11-12  hours.  Figure  42  and  Figure  43  demonstrate  the  frequency  of  different 
gap  durations.  The  vast  majority  of  gaps  in  visibility  from  Li  are  less  than  1  hr  (mainly  caused  by  eclipses 
and  the  Earth  blockage  constraint),  but  there  are  also  quite  a  few  gaps  in  the  500-600  hour  range  (caused 
by  the  phase  angle  constraint).  The  same  analysis  have  been  done  for  GEO  visibility  from  lunar  L2  (Figure 
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44-47)  and  L3  (Figure  48-51)  halo  orbit.  The  visibility  from  L2  is  almost  identical  to  Li.  The  only  difference 
between  the  Li  and  the  L2  observers  is  the  distance  and  because  L2  is  further  away  from  the  Earth  and  the 
GEO  belt  than  Li,  an  observer  at  L2  would  require  a  slightly  more  sensitive  optical  system. 

Figure  48  and  Figure  49  show  the  percentile  breakdown  of  access  and  gap  durations  from  an  L3  orbit. 
The  blue  line  marks  the  percentage  of  accesses  or  gaps  that  are  less  than  a  given  duration,  and  the  red 
line  marks  the  mean  access  or  gap  duration.  The  histogram  and  percentile  breakdown  show  the  same 
data  in  different  ways.  We  can  see  that,  for  both  accesses  and  gaps,  the  mean  value  is  not  very  useful. 
There  are  a  great  deal  of  short  gaps  and  accesses,  but  the  long  duration  gaps  and  accesses  give  us  a  more 
accurate  representation  of  a  "typical"  time  period.  The  many  short  gaps  are  caused  by  restrictions  from 
the  Earth  angle.  Roughly  28%  of  the  gaps  are  less  than  70  minutes,  the  maximum  time  a  GEO  spacecraft 
would  be  in  eclipse.  The  longer  gaps,  which  are  less  frequent,  are  caused  by  restrictions  from  the  phase 
angle.  As  seen  in  the  figure,  these  long  gaps  can  last  on  the  order  of  hundreds  of  hours.  The  average  gap 
duration  is  146  hr.  The  access  duration  has  a  similar  trend  to  gap  duration,  with  many  short  accesses  and 
fewer  long  accesses  that  can  last  for  hundreds  of  hours.  The  average  access  duration  is  42  hr. 

Another  important  consideration  for  optical  observations  is  the  distance  between  the  target  and  the 
observer.  The  distance  varies  considerably  for  an  observer  in  L3,  with  a  broad  distribution  from  1x105  km 
to  6.5x105  km.  Such  a  broad  range  of  distances  between  sensor  and  target  complicates  telescope  design, 
as  it  must  be  capable  of  varying  its  focusing  distance  by  nearly  an  order  of  magnitude  to  image  everything. 
However,  it  is  important  to  remember  that  the  longer  accesses  (which  make  up  the  vast  majority  of  total 
access  time  and  are  thus  more  desirable)  occur  when  the  distance  between  spacecraft  is  largest  (when 
the  L3  spacecraft  is  at  apogee).  Therefore,  it  may  be  possible  to  capture  most  of  the  total  access  time  by 
limiting  the  focal  distance  range  to  5x105  -  6x105  km. 

H.  Earth  Observability  Analysis  from  Libration  Point  Orbits 

Two  steps  were  defined:  Visibility  analysis  of  the  Earth  surface  using  Matlab  programs,  and  validation  of 
the  results  using  NASA's  Ephemeris  model  in  the  orbit  analysis  software  Satellite  Took  Kit  (STK).  To 
compute  cumulative  visibility,  the  Earth's  rotation  is  modeled  using  the  following  assumptions:  Earth's 
North  Pole  is  always  aligned  with  the  z-axis  in  the  CRTBP,  Earth  completes  a  rotation  in  one  sidereal  day 
(86164.09056  seconds),  and  0  deg  longitude  is  on  the  x-axis  pointing  towards  the  Moon  at  the  scenario 
epoch.  Twenty  five  hundred  locations  were  selected  on  the  halo  orbit,  each  equidistant  in  time.  At  each 
time  step,  line  of  sight  visibility  was  computed  for  all  points  on  the  latitude/longitude  grid.  Figure  52  shows 
the  average  visibility  for  two  northern  LI  halo  orbits  during  an  orbital  period.  The  most  planar  orbit,  Figure 
52(a)  provides  more  uniform  coverage  of  all  latitudes  while  the  most  vertical  orbit,  Figure  52(b),  provides 
better  coverage  of  the  North  Pole.  Figure  53  shows  the  average  visibility  of  points  on  the  Earth  for  three 
northern  L3  halo  orbits  during  one  orbital  period.  Each  halo  orbit  in  the  sequence  (a,  b,  c)  becomes 
increasingly  planar  and  more  uniform  coverage  at  all  latitudes  is  achieved.  Figures  52  and  53  show  that 
coverage  is  nearly  constant  for  a  given  latitude.  Thus,  the  next  step  in  the  analysis  is  to  compute  the 
average  visibility  for  each  latitude  for  a  variety  of  halo  orbits.  For  each  halo  orbit,  the  average  visibility  of 
each  latitude  is  plotted  with  respect  to  the  orbit's  initial  x-condition,  as  shown  in  Figure  54.  For  STK 
validation,  an  example  northern  LI  orbit  with  Az  =  45,500  km  and  an  example  southern  L3  orbit  with  Az  = 
127,800  km  were  selected.  Earth  surface  coverage  was  computed  at  a  variety  of  ephemeris  epochs.  The 
epochs  were  chosen  in  order  to  analyze  the  effects  of  the  lunar  rotational  period  (approximately  27.32 
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days)  and  the  lunar  nodal  cycle  (approximately  18.6  years).  Figures  55  and  56  show  the  average  visibility 
by  latitude  for  one  halo  orbiter  at  LI  and  L3,  respectively.  The  subplots  (a),  (b),  (c),  and  (d)  reflect  the  lunar 
nodal  cycle;  during  the  years  2015,  2020,  2025,  and  2029,  the  Moon's  right  ascension  of  the  ascending 
node  with  respect  to  the  Earth  is  approximately  180,  90,  0,  and  270  deg,  respectively.  For  each 
year, sixteen  epochs  were  selected.  Each  subsequent  epoch  was  6.83  days  (one-fourth  of  the  lunar 
rotational  period)  after  the  previous.  Thus,  Epoch  2  is  6.83  days  after  Epoch  1  and  so  forth.  For  reference, 
the  nominal  coverage  computed  in  the  CRTBP  is  also  shown  in  each  subplot.  Also  note  that  the  orbit's 
initial  state  in  the  CRTBP  is  on  x-axis  with  positive  y-velocity.  Figures  55  and  56  show  that  the  average 
coverage  by  latitude  approximately  repeats  every  fourth  epoch,  correlating  to  the  lunar  rotational  period. 
The  variations  between  every  fourth  epoch  are  slightly  larger  for  LI  orbits  than  for  orbits  at  L3. 

Also,  coverage  is  computed  for  single  satellites  and  constellations  of  spacecraft  in  halo  orbits  about  LI  and 
L3.  The  goal  is  to  achieve  the  highest  surface  coverage  using  a  minimum  number  of  spacecraft. 
Constellations  considered  include  a  two-spacecraft  constellation  with  one  LI  orbiter  and  one  L3  orbiter, 
a  three-spacecraft  constellation  with  one  LI  orbiter  and  two  L3  orbiters  in  mirrored  northern  and  southern 
orbits,  and  a  four-spacecraft  constellation  with  two  LI  orbiters  and  two  L3  orbiters.  Figures  57  and  58 
show  the  average  percent  visibility  of  the  Earth's  surface  by  at  least  one  spacecraft  and  average  percent 
visibility  at  each  latitude  for  the  constellations  that  provide  the  best  and  worst  coverage.  For  STK 
validation,  the  orbits  that  provided  the  best  coverage  were  analyzed  in  STK  for  the  epochs  used  in  the 
single-spacecraft  coverage  analysis.  For  the  three-spacecraft  constellation  the  coverage  from  the  full 
ephemeris  is  compated  with  that  from  the  Matlab  analysis  in  Fig.  59.  Similar  to  the  coverage  provided  by 
one  orbiter.  Fig.  59  shows  that  the  average  coverage  by  latitude  approximately  repeats  every  fourth  epoch 
for  the  constellation  while  the  variations  between  different  epochs  are  on  a  much  smaller  scale.  It  is  also 
observed  that  the  coverage  in  the  ephemeris  is  better  at  the  poles  than  in  the  Matlab  analysis  of  the 
circular  restricted  three-body  problem.  Still,  the  CR3BP  is  a  reasonably  good  model  for  the  coverage  that 
can  be  expected  in  the  ephemeris.  Furthermore,  the  fact  that  certain  L3  southern  orbits  may  provide 
more  coverage  of  the  Northern  hemisphere  is  nullified  by  the  fact  that  two  reciprocal  orbits  are  used  at 
L3  for  the  constellation. 
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Figure  3-Several  example  orbits  from  the  family  of  halo  orbits  about  the  Earth-Moon  L3  point. 


Figure  4-Values  of  i  around  a  periodic  orbit. 
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Figure  5-Stable  (green)  and  unstable  (red)  exterior  manifolds  of  an  LI  halo  orbit. 
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Figure  6-Stable  manifold  trajectories  of  two  LI  LPOs  with  significantly  different  stability  coefficients 
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Figure  7-The  perigees  of  WS  trajectories  of  an  L3  halo  orbit  with  C=2.42  propagated  for  one  year  (left)  and  two  years  (right). 
The  color  of  the  marker  indicates  the  day  on  which  perigee  occurs. 


Figure  8-Planar  views  of  the  stable  (left)  and  pseudo-stable  (right)  manifold  trajectories  which  depart  from  the  same  location 

on  the  nominal  L3  halo  orbit. 
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Figure  9-Optimal  two-impulse  transfer  from  LEO  to  a  southern  LI  halo  orbit  with  C=3.10.  Delta-v2  is  indicated  by  the  pink 
circle.  The  genetic  algorithm  scatter  plot  is  shown  converging  to  the  optimal  solution 
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Figure  10-Optimal  two-impulse  transfer  from  GEO  to  a  southern  LI  halo  orbit  with  C=3.10.  Delta-v2  is  indicated  by  the  pink 
circle.  The  genetic  algorithm  scatter  plot  is  shown  converging  to  the  optimal  solution. 


Figure  11-Direct  and  an  alternative  locally  optimal  transfer  from  LEO  to  a  southern  LI  halo  orbit  with  C=3.10.  Delta-v2  is 

indicated  by  the  pink  circle. 
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Figure  12-Direct  and  an  alternative  locally  optimal  transfer  from  GEO  to  a  southern  LI  halo  orbit  with  C=3.10.  Delta-v2  is 

indicated  by  the  pink  circle. 
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Figure  13-Effects  of  thrust  magnitude  and  direction  dispersions  at  manifold  injection  on  miss  distance  and  transfer  times 
without  guidance  for  LEO  transfers.  Miss  distance  dispersions  are  for  thrust  errors  in  magnitude,  in-plane  direction,  out-of- 

plane  direction 
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Figure  14-Effects  of  thrust  magnitude  and  direction  dispersions  at  manifold  injection  on  miss  distance  and  transfer  times 
without  guidance  for  GEO  transfers.  Miss  distance  dispersions  are  for  thrust  errors  in  magnitude,  in-plane  direction,  out-of¬ 
plane  direction 
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Figure  15-  LQR-based  guidance  with  manifold  injection  errors  for  LI  halo  orbit  transfer  from  LEO. 
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Figure  16-LQR-based  guidance  with  manifold  injection  errors  for  LI  halo  orbit  transfer  from  GEO. 
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Figure  17-LQR-based  limited  thrust  guidance  (Tm  =  0.01  mm/sA2)  with  manifold  injection  errors  for  LI  halo  orbit  transfer 

from  GEO. 
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Figure  18-Liapunov  exponents  for  the  LI  halo  orbit  stable  manifold  associated  with  the  optimal  transfers  from  LEO  and  GEO 

in  Figs.  9-10. 
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Figure  20-lmpulsive  guidance  with  manifold  injection  and  TCM  thrust  errors  using  GEO  departure. 
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Figure  21-Direct  and  spiral  departure  transfer  trajectories  from  GEO  to  a  LI  Liapunov  orbit  with  corresponding  continuous 
thrust  magnitudes  using  mixed  impulsive  and  continuous  thrust  and  quadratic  performance  index  in  the  OCP. 


Figure  22-Direct  and  spiral  departure  transfer  trajectories  from  GEO  to  a  LI  Liapunov  orbit  with  corresponding  continuous 
thrust  magnitudes  using  mixed  impulsive  and  continuous  thrust  and  minimum  fuel  performance  index  in  the  OCP. 
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Figure  23-A  pseudo-manifold-based  transfer  trajectory  from  a  185  km  LEO  to  an  L3  halo  orbit  with  C  =  2.40207. 
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Figure  24-Perigee  altitudes  (km)  for  different  trajectories  with  respect  to  the  insertion  point  into  the  L3  halo  orbit,  where 

delta-v2=100  m/s. 
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Figure  25-Different  families  of  transfers  from  a  185  km  altitude  LEO  to  the  L3  halo  orbit. 
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Figure  27-Transfer  trajectory  to  the  L3  halo  orbit  which  has  the  smallest  possible  insertion  maneuver.  C3=-2.05  kmA2/sA2, 

AV2=0.015  km/s,  TOF=119.6  days. 
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Figure  28-Transfer  trajectory  to  the  L3  halo  orbit  which  has  the  least  possible  time  of  flight.  C3=-2.02  kmA2/sA2,  AV2=0.09 

km/s,  TOF=31  days. 
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Figure  29-Visibility  of  Earth  surface  using  a  spacecraft  on  LI  halo  orbit 
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Figure  30-Visibility  of  geosynchronous  orbit  using  a  spacecraft  on  LI  halo  orbit 
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Figure  31-Average  percentage  of  Earth  surface  visibility  by  one  or  more  satellites  in  a  constellation  of  two  satellites  in  L3  halo 

orbit  and  one  on  LI  halo  orbit 
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Figure  32-Screenshot  of  STK  showing  the  LI,  L2,  and  L3  libration  point  orbits  used  in  the  analyses  presented  here  and 

illustrated  in  the  Earth-Moon  rotating  coordinate  frame 


Figure  33-Portion  of  time  where  target  is  visible  for  various  altitudes  and  inclinations 
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Figure  34-Mean  access  duration  for  various  altitudes  and  inclinations. 
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Figure  35-Mean  gap  duration  for  various  altitudes  and  inclinations. 
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Figure  36-Effect  of  each  constraint  on  percent  time  visible  for  the  LEO  and  MEO  regimes. 
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Figure  37-Percentage  of  GEO  spacecraft  visible  from  LI  over  time 
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Figure  38-Percentage  of  GEO  spacecraft  visible  from  L2  over  time 
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Figure  39-Percentage  of  GEO  spacecraft  visible  from  L3  over  time 
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Figure  40-Breakdown  of  access  durations  from  LI. 
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Figure  41-Histogram  of  access  durations  from  LI 


Figure  42-Breakdown  of  gap  durations  from  LI 
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Figure  43-Histogram  of  gap  durations  from  LI 


Figure  44-Breakdown  of  access  durations  from  L2. 
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Figure  45-Histogram  of  access  durations  from  L2 


Figure  46-  Breakdown  of  gap  durations  from  L2 


Figure  47-Histogram  of  gap  durations  from  L2 


Figure  48-Breakdown  of  access  duration  from  L3. 
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Figure  49-Histogram  of  access  durations  from  L3 


Figure  50-Breakdown  of  gap  duration  from  L3. 
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Figure  51-Histogram  of  gap  durations  from  L3. 


(a)  x0  =  316.839  km.  C  -  3.17  (b)  x0  -  359.162  km.  C  =  2.98 


Figure  52-Percent  visibility  for  latitude/longitude  points  on  the  Earth  surface  for  selected  northern  LI  halo  orbits  assuming 

Earth's  spin  axis  is  identified  with  z-axis  of  CR3BP. 
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Figure  53-Percent  visibility  for  latitude/longitude  points  on  the  Earth  surface  for  selected  northern  L3  halo  orbits  assuming 

Earth's  spin  axis  is  identified  with  z-axis  of  CR3BP. 
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Figure  54-Average  Earth  visibility  by  latitude  for  northern  LI  and  L3  halo  orbits. 


(b)  2020 


(c[  2025 


(d)  2020 


Figure  55-Average  visibility  by  latitude  at  different  ephemeris  epochs  for  an  LI  orbiter  from  STK  analysis. 


Figure  56-Average  visibility  by  latitude  at  different  ephemeris  epochs  for  an  L3  orbiter  from  STK  analysis. 
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(a)  Two spacccraA  corateDatum:  Northern  L.j  and  (h}  Three-  apweetraft  conatdlation:  Northern  Lj  . 
southern  Lg  northern  Lj,  and  sautheni  L3 


Figure  57-Average  percent  visibility  of  the  Earth's  surface  by  at  least  one  spacecraft  in  a  constellation  of  a)  two  spacecraft 

and  b)  three  spacecraft  in  LI  and  L3  LPOs. 
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Figure  58-Average  percent  visibility  at  each  latitude  for  the  constellations  that  provide  the  best  and  worst  coverage  in  a 

constellation  of  a)  two  spacecraft  and  b)  three  spacecraft. 


Figure  59-Visibility  by  latitudes  at  different  ephemeris  epochs  for  a  three-spacecraft  constellation:  comparison  of  STK 

analysis  (Epochs  1-16)  with  the  nominal  CR3BP  analysis. 


